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Abstract 

In this paper, we study a curvaton model where the curvaton is acted by Galileon field. We cal- 
^v^j ' culate the power spectrum of fluctuation of G-curvaton during inflation and discuss how it converts 

to the curvature perturbation after the end of inflation. We estimate the bispectrum of curvature 
perturbation induced, and show the dependence of non-Gaussianity on the parameters of model. It 
j ^ ■ is found that our model can have sizable local and equilateral non-Gaussianities to up to O(10 2 ), 

which is illustrated by an explicit example. 

*~ 1 . PACS numbers: 

I. INTRODUCTION 

Making our universe expand fast enough, it can naturally solve many notorious problems brought out by hot Big Bang, 
such as flatness problem, horizon problem, monopole problem and so on. Moreover, during inflation the quantum 
fluctuations generated at the initial stage can be stretched out of the horizon to form classical perturbations, which 
can provide seeds for the formation of the structure of our universe. An inflation model also succeeds in producing 
nearly scale-invariant power spectrum of scalar perturbation and tiny gravitational waves, which fits very well with 
today's observational data The non-Gaussian corrections of perturbations during inflation can also be large or 
not, according to various inflation models, which is waiting for constraints from new and more accurate data in the 
near future @. 

Usually, the perturbation of inflation is generated by the inflaton field itself, which is the simplest way to have 
curvature perturbation. However, it is not the only choice. Perturbations generated in such a way depends on the 
potential of the inflaton field, and thus puts very severe constraint on inflation models. In order to relax such a 
constraint, as Lyth and Wands have pointed out, perturbation can also be generated from another field that has 
nothing to do with the inflaton field, namely, the curvaton Q, see also relevant works on curvaton mechanism in 
Refs. @, § and earlier [t| [l(| ■ Curvaton field is usually assumed to be a scalar field with light mass and decoupled 
from all the other kinds of perturbations, thus the perturbation produced by curvaton can be independent on the 
nature of inflation. Moreover, since the curvaton is subdominant during inflation, it can only produce isocurvature 
perturbation. This isocurvature perturbation has to be converted into curvature perturbation at the end of inflation, 
so it depends on what happened after the inflation terminated. Usually, there are two cases in which this conversion 
can be available: First, when the inflaton decays into radiation after inflation, the curvaton field becomes dominate; 
Second, the curvaton field decays as well before its domination, and reaches equilibrium with radiation that decays from 
inflaton. According to different case, the amount of the curvature perturbation converted from curvaton perturbation 
may be different. Curvaton scenarios have been widely studied in, for instance, 

In the original curvaton paper, it was suggested that curvaton is made of a canonical scalar field. However, other field 
models can be considered to act as a curvaton. People have considered curvaton of Pseudo-Nambu-Goldstone-Boson 
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Pj , DBI-type flp o r its curvaton brane implement [14], multi- field [Hj], lagrangian multiplier field [l6j |. Horava- 
Lifshitz scenario [13] and so on. Very recently, a new kind of models has been proposed and studied extensively, 
which is called "Galileon" models id- The original version of these models is a generalization of an effective field 
description of the DGP model [l9| . These models includes high derivative operator of the scalar field, however, due 
to some "delicated design", its equation of motion remains second order, 1 thus it can violate the NEC without 
incorporating any instability modes. Due to such an interesting property, such kind of model jrenerically can admit 
superluminal propagation [2lfj and perhaps even closed timelike curves (CTCs) [22j (see also [23| for CTCs without 
Galileon) . Moreover, these models can be applied onto various evolution period of our universe, such as dark en ergy 
24, HH, inflation [26l - [28j j . reheating [2{|, bouncing [3(J, the slow expansion scenario f3ll o f primordial universe [32j . 
33 1, and so on. As an extension, these models can also be generalized to DBI version [3J], the K-Mouflage scenario 
35 j, the supersymmetric Galileon (36|, the Kinetic Gravity Braiding models [24]], the generic Galileon- like action 
20t l37l l3c| and others. See also e.g. Refs. (39l-l46| for various study of their phenomenologies. 

In the present work, we study the scenario where a Galileon field behaves as a curvaton, which we dubbed as 
" G- Curvaton!^ scenario. Due to the higher derivative term, Galileon is expected to have some features on generating 
perturbations, such as getting large tensor-scalar ratio in "G-inflation" scenario. Our paper is organized as follows: in 
Sec. II we review the original curvaton scenario, by taking curvaton to be the simplest one, i.e., the canonical scalar 
field. In Sec. Ill we study our "G-Cuvaton" model. We first investigate its perturbation, obtaining a scale-invariant 
power spectrum. Then we discuss how it converts to the curvature perturbation at the end of inflation, considering 
both cases where curvaton decays before and after it dominates the universe. We also study the non-Gaussianitics 
generated from our model, both local type and equilateral type. Finally we present an explicit example to show how 
the observable quantities could be effected. Sec. IV comes our conclusion and discussions. 



II. REVIEW OF THE SIMPLEST CURVATON MODEL 



In this section, we would like to review how the mechanism works for the simplest curvaton model, which is made 
of a canonical scalar field. The Lagrangian for the curvaton field is: 



C tT = -~d fl ad»a-V(a) . (1) 



Note that here we are using the metric with notation ds 2 = —dt 2 + a 2 (t)Sijdx l dx : ' . As for curvaton, the effective mass 
of a needs to be very light, which put the constraint \V acr \ <C H 2 to the potential V(a), where V a(T = d 2 V/da 2 and 
H denotes the Hubble parameter of the universe. Moreover, during inflation H is almost a constant. Thus in general, 
one may define the parameters 

H 2 3H 2 ' ( ) 

which are required to be far smaller than 1 and slowly varying. Moreover, since curvaton is subdominant part of 
the universe during inflation, the energy scale of the potential should also be lower than that of inflation, namely 
V(a) « 3M» H\ 

Suppose the field generates fluctuation during inflation, namely a{x) — > cr (t) + (5cr(i,x), then from Eq. ([T]) we can 
easily get the equation of motion for the field fluctuation 5a as: 

5a k + 3H5a k + {{k/af + V aa )5a k = , (3) 

where 5a k is the Fourier presentation of 5a with momentum mode k, and a is the scale factor. Note that here we have 
already neglected the coupling of 5a k to the metric perturbation, as has been done in [6] . Thus the power spectrum 
of the field fluctuation 5a k at the horizon crossing is given by 

'•-sfI 4 *!' -(£)'■ (4) 

where the star denotes the time of horizon exit, k = a*iJ*. The spectral index of the spectrum can also be given by: 

n.-l=^ = -2e + 2e«l, (5) 



1 This idea is actually pioneered by Horndeski thirty years ago, in [20| . 
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which shows the nearly scale-invariance of the spectrum V a - 

However, since the curvaton is the subdominant part of the universe during inflation, the perturbation generated 
by it can only be of isocurvature type. This type of perturbation can only be converted into curvature perturbation 
when curvaton dominates, or become equilibrium with other part of the universe. Either these cases will not happen 
during inflation, however after inflation, when inflaton field decays into radiation whose energy density may decrease 
more rapidly than curvaton field, both the two cases will happen, depending on when the curvaton field will decay. 
In the case when the curvaton decays late, it will exceed over the radiation decayed from inflaton and dominate the 
universe, however in the case when the curvaton decays early, it may become equilibrium with the radiation. The 
curvature perturbations for a given matter with energy density p in spatial- flat slicing is given by [47] : 

C = -H 6 4 , (6) 
P 

and the separately conserved curvature perturbations for radiation and curvaton field therefore read: 

Cr = -H 8 ^ = ^ , (7) 
Pr 4 p r 

<« = -H 6 -^= w ^- ] , (8) 

Pa 3{Pa+Pa) 

respectively, where p a is the pressure of the curvaton field. Using these results, and assuming that the isocurvature 
perturbation convert to curvature perturbation instantly, then the curvature perturbation generated by such conversion 
reads: 

£ = 4p r (r + 3(p CT + Pa)(g ^ 3(Pa + Pa)(a , g , 
Ap r + 3(p a + Pa) ^Pr + i{Pa+Pa) 

where in the last step we neglected the curvature perturbation of radiation, £ r . Define the energy density ratio of a 
over radiation, r = p a / p r , then in the first case where the curvaton dominates, we have r> 1, Eq. (|9]) becomes: 

C - Cr , (10) 

while in the second case where a become equilibrium with radiation, we have p a <C p r , Eq. ([9]) becomes: 

C ~^r(l + OG , (11) 

where w a = Pa/pa is the equation of state of the curvaton field. 

Furthermore, we investigate the non-Gaussianities of the perturbation generated by the curvaton field. The local 
type non-Gaussianities of curvature perturbation are given by: 

C = C 9 + \f£l al C 2 9 , (12) 

where the subscript "g" denotes the Gaussian part of £ while Jml is the so-called nonlinear estimator. For local type, 
/yv£ a ' can t» e estimated by using the so-called 5N formalism, where SN is the variation of the number of efolds N of 
inflation [Hj]: 

C = SN = N, a 5a + ^N aa Sa 2 + ... , (13) 
where A r cr...cr = d n N/da n and the same notations hereafter. Comparing Eqs. (fT2"|) and (TIB")) one can easily find that 



In c l 1 can t> e presented using N a and N namely, 



local 
NL 



^ (14) 
6 Nl [ ' 



For nonlocal type, however, things will become a little bit more complicated, since non-Gaussianities also exists 
in the field fluctuation 5a itself, and SN formalism will not be valid any longer. In this case, we could express the 
3-point correlation function of £ as: 

(|C(*i)C(*a)C(*3)|) = (27r) 3 <5 3 (^ k i )B(k 1 ,k 2 , fca) (15) 
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where B{k\, &2, ^3) is the shape of the non-Gaussianties and 



(|C(fci)C(fc 2 )C(fc 3 )l> = -*T / £ft , <|[C(*,*i)C(* J *2)C(t,*3),^S.t(* , )]l) . ( 16 ) 

Jt 



with "Hf nt being interaction Hamiltonian of curvaton in momentum space. The nonlinear estimator is defined as: 

10 IlLikf B{kxMM 



£ nonlocal 
JNL 



(17) 



C 3 (2tt) 4 £•=! fcf [PdWPdh) + 2perm.} ' 

III. THE G-CURVATON MODEL 

In this section, we study our model of which the curvaton field is of Galileon type. For simplicity but without 
losi ng g enerality, we consider the Lagrangian of curvaton has the general third order Galileon field, which is given by 

Mm 

Sac = J d 4 x^g~[K(a,X) - G(a,Xpa] , (18) 

where K and G are generic functions of a and X = — d^ad^a/2 is the kinetic term of the field a. The Lagrangian 
used is inspired by the original Galileon construction but does not respect the Galileon symmetry, which was called 
'Generalised Galileons' by Deffayet et al.[.24j. Note that more generalized Galileon model containing higher order 
operators of or the couplings of a to the gravitational part was constructed in Ref. [37]. From action (fT5|) . the 
energy-momentum tensor T^ v has the form of: 

Tfu, = K^d^ad^a + Kg^ u - d^Gd^a - d^Gd^a + g flu d x Gd x a - G^od^ad v a . (19) 

Taking the homogeneous and isotropic background where T^ v has the form of diag{p ai a 2 (t)p tT ,a 2 (t)p a ,a 2 (t)p a -}, we 
can furtherly obtain the energy density and pressure of a as: 

Pa = 2K, X X -K + 3HG, x <j 3 - 2G, a X , (20) 
Pa = K -2{G, a + G. x a)X . (21) 

Moreover, by varying action (|18[) with respect to a, we obtain the equation of motion for a: 

K tX Ua - K iXX (d„d v a)(d fi ad u a) - 2K t x<rX + K, a -2 (G t<7 - G, X *X) Ua 
+G.x [(c^cUXc^V) - (Ua) 2 + R^aff'a] + 2G,xa(d^a)(d^ad"a) 
+2G t<7t7 X - G tX x (d»d x o - g^Ua) (d^a) d v ad x a = . (22) 

which can be simplified in FRW universe where the background are homogeneous and isotropic: 

K x (a + 3Ha) + 2KxxXa + 2K X *X - K, a - 2 (G CT - G, X *X) (a + 3H&) 

+6G, X [(HXy + 3H 2 X] - AG^xaXa - 2G^ <J X + QHG tX xXX = . (23) 

A. Scale-Invariant Power Spectrum for Field Fluctuation 

To study the perturbations of the curvaton field, we can split the scalar field a into a[x) — > <Jq{£) + <5ct(£,x), 
where ao represents the spatially homogeneous background field, and the 8a stands for the linear fluctuation which 
corresponds to the isocurvature perturbation during inflation. Taking the spatial-flat gauge and using the assumption 
that there is no coupling between 5a and other perturbations, one can have the equation of motion for field fluctuation 
in momentum space Sa^ as: 

Sa k + f 3 + HSa k + ^-5a k + M 2 eff 5a k = , (24) 
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where we have defined c 2 = C/T>, with 



C = K tX + 2G,x((To + 2Ha ) + 2G,xxXao-2(G t(T -G,x<rX) , (25) 
V = K. x + 2K. XX X + &HG >X & - 2 {G >a + G. Xa X) + 6HG, xx Xa Q . (26) 

Note that in order to avoid ghost or gradiant instabilities in our model, one must require C > 0, T> > 0, which leads 
to the non-negativity of the sound speed squared: c 2 > 0. The effective mass squared is: 

M\f f = ^ t [a 3 {K, XlT & + 3HG, x <rCr% - G, a(7 a )] - K, a(T + G, aa Ua . (27) 
For later convenience, we also introduce the following "slow variation" parameters: 

H d s V M^f 

^ - - W , *g - , S G — , & - -g^-, (28) 

which are assumed to be small but not neglected. When reduced to the simplest curvaton model given by ([1]), eg = e, 
= f and s G = S G = 0. 

For solving Eq. (|24[) and computing the power spectrum, it is convenient to turn to the conformal coordinate where 
conformal time r is defined as dr = dt/a. Using a new variable Uk = zSa^ where z = a^/V, the equation of motion 
can be written in the Fourier space as 

< + (elk 2 -^j Uk = 0, (29) 

and the prime denotes differentiation with respect to r. Under the slow roll approximation where all the parameters 
in (|28|) are small, we find: 

z" 2 

Z T z 

Thus we finally can obtain the power spectrum of Sa-^ at horizon crossing as (using definition (|3])): 

H 2 

Vs °* = 4^ (31) 

with the spectral index given by 

d In *P 

n a -\ = — - -2e G - 3s G - S G + 2& « 1. (32) 
dmk 

From Eq. (|32[) we can see that the isocurvature perturbation generated by our curvaton field can give rise to a very 
flat power spectrum which is nearly scale-invariant. Moreover, comparing to usual curvaton case the amplitude 
of the power spectrum is suppressed by Galileon-like nonlinear terms such as T>, however, when G(cr, X) = and 
K(a,X) — X — V(a), the result exactly reduces to that of usual curvaton case. When G(a,X) = and K{a, X) is 
DBI-typc, the result is that of DBI-curvaton[l3j]. The result given here is also consistent with the case where Galileon 
field act as an inflaton field and generate curvature perturbations [26j] 2 , and the similar property has also been found 
in other featured inflation models, such as DBI inflation [491 . 



B. Generating the Curvature Perturbation from Curvaton 



From last paragraph we learned that our G-Curvaton model is able to give rise to perturbations with nearly scale- 
invariant power spectrum, however these perturbations are of isocurvature ones. As has been mentioned before, 
curvature perturbations can be obtained after inflation, so now we consider the epoch when inflation has ceased and 
the inflaton has already decayed to radiation, during which the universe is filled with the curvaton field p a and the 



2 Strictly speaking, as authors of Ref. pointed out themselves, what they calculated is not usual comoving curvature perturbation 
but another variable which coincide with comoving curvature perturbation only in large scales. 
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radiation p r . From this moment on, the isocurvature perturbations began to convert into curvature ones, and this 
conversion will complete in two possible cases, namely, either when the curvaton dominates the universe, or decay 
as well, whichever is earlier 3 . We assume that the conversion as well as curvaton decay happens instantaneously, so 
that we can separately consider the curvature perturbations for each component (radiation and curvaton) and just 
use Eq. © to calculate the final total curvature perturbation of the universe. From Eqs. flSJ) and ([9]) we know that 
the final curvature perturbation can be expressed as: 

C= , - ; , (33) 

where the density perturbation Sp a can furtherly be expanded with respect to 8a. The linear order of 5p a is given by: 

(5 (1) p<r ^ p<j,<j5(t 

= {2K >Xa X - K >!T + 3HG, Xa a% - 2G, aa X)6a , (34) 

while the second order of 5p a reads: 



5 (2) Pa 2t -Pcp^G 1 



1 

2 1 

= i {2K >Xaa X - K, aa + 3HG, Xaa a 3 ~ 2G^ aa X)5a 2 (35) 

respectively. 

Now we can consider the two cases separately. First, if the curvaton dominates the energy density before decays, 
we can use Eqs. {§]), <TlO]) as well as (j34l) to have the final curvature perturbation as: 

3(pa+Pcr) 3(pcr+Pcr) 

1 / K Xa &% - K ia + ZHG, Xa dl - G, aa a 2 
3 ^ K x a 2 + 3HG <x a$ - 2G. a a 2 _ G , x a Q a 2 



6a , (36) 



and the power spectrum of curvature perturbation is: 



V (J) = Jf_| C |2 

C ~ tor*™ 

■2 



1 / H 2 \ ( K, Xa a 2 - K, a + ZHG, Xa dl - G^a 2 
9 \^alv) \K X &$ + ZHG tX al - 2G> 2 - G :X a a 2 



(37) 

where we also made use of our previous results pip . Second, if the curvaton decays before its dominance, it will only 
contribute part of the energy density of the universe with some fraction r as defined before. Using Eq. in this 

case the curvature perturbation is given by 

An) _ r_5^_pa_ _ rp ff|fffo 



r 



\Xaa Q — JX tlT -f on<^ iXtT a% — G. aa a 



4\ K iX al-K + 2>HG tX al-G ia &l 
and the curvature perturbation power spectrum is: 



5a , (38) 



v (ii) = r 2 ( K x „&* - K g + 3HG, XtT a 3 - G ^a 2 \ 2 fc 3 2 
< 16 V, K iX a 2 -K + 3HG iX a 3 -G ir7 a 2 J 2^ 1 1 

( H 2 \ ( K, Xa a 2 - K t<T + 3HG, x<T a 3 - G, aa a 2 
16 \4w 2 c^Vj \ K, x a 2 -K + 3HG, x a$-G, a a 2 



(39) 



Actually, besides the standard model radiation, the curvaton can also decay into other products such as dark radiation or dark matter. 
In that case, there may have residual isocurvature perturbations and the non-Gaussianities may also be different |50( (see also |5l| for 
a review). In this paper we assume that curvaton decays to standard model radiation for simplicity. We thank A. Mazumdar for point 
out this for us via private communication. 
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Here the supscripts (7) and (77) refer to the two cases respectively. 

From Eqs. (j3"7| and (f3T?)) one can see that comparing to usual curvaton case, the amplitude of curvature perturbation 
will also be modified with a different prefactor, which contains contribution from the high derivative term G(a, X). 
So one can naturally expect G-Curvaton model has some feature that could be distinguished from usual K-essence 
curvaton models. 



C. Non-Gaussianities generated by G-Curvaton 



In this paragraph, we extend our study to non-linear perturbations of our model, i.e., non-Gaussianities. First 
of all, we consider the local type non-Gaussianities generated by the curvaton model. For local type, the nonlinear 
parameter f x c L al has been defined in (Tl2|) . From Eqs. (|33H35p and using the SN formalism (|13|) . we can also express 



yioca; m terms of energy density as: 



elocal 
INL 



C 



~ [Ap r + 3(p a + p a )} 

6 Per. (7 



(40) 



Now we consider the two cases separately. For the first case where curvaton dominates the energy density before 
decays, one gets: 



ocaL 
NL 



(I) 

c 



5 {Per + Pa)Po 
2 



(K x & 2 + 3HG tX a 3 - 2G> 2 - G iX a a 2 )(K tXaa & 2 - K. aa + 3HG, x<T<T a 3 - G^ a a 2 ) 



(K, Xa a 2 - K. a + 377G Xtr d% - G^o 2 ) 2 



(41) 



and for the second case where the curvaton decays and never dominates the energy density, we have 

(m _ io PaPa , aa 

c 



NL 



3 r Pi, 



10 
3^ 



{K. x a 2 -K + 3HG, X & 3 - G ^b 2 )(K. Xaa b 2 - K. aa + 3HG, Xaa & 3 - G^a 2 ) 



(K, Xt7 a 2 - K M + 3HG^ Xa & 3 - G, aa <r 2 ) 2 



(42) 



respectively. 

Then we turn to the nonlocal type of non-Gaussianities. The nonlocal type of non-Gaussianities is more complicated, 
as described in Eqs. (|15H17[) . The interaction hamiltonian is H^ nt = —£3, where £3 is the 3-rd order perturbed 
Lagrangian. Starting from action (|18p and after straightforward but rather tedious calculation, we get the interaction 
hamiltonian at the leading order with respect to slow roll parameters: 



77 p ~) 



d 3 k 1 d?k 2 d 3 k 3 , 



(2 7 r) 3 ,5 3 (k 1 + k 2 + k 3 ){+ai!(k 2 • k 3 )Sa(t, ki)*<r(t, k 2 )Sa(t, k 3 ) + a 3 L 2 Sa(t, ki)Sa(t, k 2 )5a(t, k 3 



+a7 3 (k 2 • k 3 )5a(t, kx)Sa(t, k 2 )Sa(t, k 3 ) + a~ l L A (ki • k 2 )k 2 3 Sa{t, kx)5a(t, k 2 )Sa(t, k 3 )} , 
where we have defined 



(43) 



Li 
L 2 

L 3 

Li 



G 



-(5HG yXX0 &l + 277(7 xo + -^K.xxxoO'o + HG yXXX0 &Q + K yXXO a Q ) 
— -^{K,xxo&o + 2G, xxo & <jo + 377G i xxoo'o + 477GX0) , 

~l G ' X0 ■ 



(44) 
(45) 

(46) 

(47) 



Moreover, from the solution of the canonical variable Uk we can obtain the mode solution of the perturbation 
variable (in conformal time) S<j(t, k) as: 



Sct(t, k) = q(r, k)a k + q*(r, k)al k , q(r, k) 



111 



: (1 + ic s kr)e 



-ic s kr 



(48) 
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where and al k are production and annihilation operators respectively. By substituting Eq. (|43H48[) into (1161) , one 
can get the result of 3-point correlation function (\8cr{ki)8a(k2)5a(k 3 )\) . Note that when carrying out integrations 
with respect to conformal time r, we assume that the Li's (i — 1, 2, 3, 4) are all slow varying and can be roughly taken 
out of the integrations, which greatly simplifies our calculation. Similar assumptions has been made when calculating 
more general single field inflation with nonminimal coupling to Gravity as well as Gauss-Bonnet terms [43[ though 
more rigid calculation with nonminimal coupling only was performed in |52j . The final shape of (\5cr(ki)5a{k2)5a{hi)\) 
than reads: 



Bsa(h, fa, fa) 



3Li# 6 



3LoH 5 



s 1 Z J i vjtj ift 

(2 k 2 k) - 7 k\k) + 4 ]T kik$ - lSkfkikl £ fa 



V 3 c 8 s fafak 3 K 3 V 3 c 6 s fafak 3 



i^j i^j i& 

-Afafafa J2 k \fa ~ 2Afafafa J2 kitf + \2fafafa Y, k t + Yl k ^ 



(49) 



and the corresponding non- linear estimator f^L can be calculated according to Eq. (|17[) . As an example, here we 
consider the equilateral case where fa — fa = fa, in which the estimator then reduces to: 



equil 
NL 



L 2 H 17L 3 H 13L 4 H 2 



8a 27"9Vc 2 s 92? 3625c;? l8Vc 4 s ' ' 
while using SN formalism (|13p . the nonlinear estimator for curvature perturbation reads: 



equil 
NL 



C 



= N, a f 



equil 
NL 



From Eqs. (|3"3")l and ([M]), we can easily get N_ a = p G ^l\Ap T 



5(7 



(50) 



(51) 



which gives 



equil 
NL 



P<7,<7 requil 

JNL 



3(p CT +p a 



L X H 2 L 2 H 17L 3 H 13L 4 H 2 



81(p a +p a ) K 9Vc 2 92? 
for the case of which curvaton dominates before decays, and 

(ii) 



mVcl 18X>cf 



(52) 



requil 
JNL 



'P(T,<7 requil 

~ ~ipV fNL 

hrp a „ ,LiH 2 L 2 H 17 L 3 H 13L A H 2 



54p a v 9Vc 2 s 9V 
for the case of which curvaton decays and never dominates. 



3625c? 



18X>cf 



(53) 



D. A concrete G-Curvaton model 



In previous parts of this section, we have presented the whole process of how our G-Curvaton model works, including 
its background evolution, scale-invariant isocurvature perturbation generation, curvature perturbation conversion as 
well as different types of non-Gaussianities. However, due to the involvement of our model, especially containing 
high derivative terms, the above general analysis can only be rather qualitative. Moreover, unlike the usual curvaton 
model, there are many uncertainties in our model with general form (fTHJ) and it can have various decaying mechanisms, 
each of which may have different results for curvature perturbations and non-Gaussianities. In order to make things 
more specific, it is necessary to focus on some explicit models to see how our models can be different from the usual 
curvaton models. 

Before constructing models, let's investigate how many constraints we have to consider for our model. First 
of all, as was mentioned in previous section, a curvaton model must have light mass, that is, -M 2 ff *C H 2 . In 
slow roll approximation, we have a <C Afi, so the first term of Eq. (|27|) can then be neglected, which makes 
Mljf ~ — K i<Ta + 3HG, aa &. Therefore, if we choose \K i<TIT \ and \G iCrcr \ to be small enough, it will be safe for curvaton. 
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Another way is both K >aa and 3H&G. aa may be large in amplitude, but are of similar value and opposite sign. In this 
case, they can be canceled to have a relatively small value, which may need some fine tuning in the model. The second 
constraint comes from the observations. The current observation data gives very tight constraint on the amplitude 
of the curvature perturbations, for example, the WMAP-7 measurement of the CMB quadrupole anisotropy requires 
Vq ~ 2.4 x 10~ 9 [5j]. In order for the amplitude of the curvature perturbation in (|37j) and (f39|) to meet the data, one 
can furtherly constrain the form of Lagrangian of our model, namely K{a, X) and G(a, X) and their field dependence. 

Taking account to both constraints from above, we can consider that our model may have the Lagrangian with 
form of: 

K(a, X) = X - V(a) , G(a, X) = -g(a)X (54) 

as an example. Here we require both V aa and g^ aa be small enough to give rise to small effective mass needed for 
curvaton. For background evolution, from eqs. ((20]), (|2"Tj) and (|23|) . we have the following equations: 



PcT = X(l - 6gH& + <?>o) + V , (55) 

Pa + Pa = 2X(1 + ga Q - 3gH&o + g, a &l) , (56) 

cto + 3H&o + 2g t<T &i(To + \g^aa^t - 3gH&Q - 6gH& ^o - §gH 2 al + V a = . (57) 
For perturbations, from Eqs. (|2"5"|) and ([2"Bf we can get: 

C = 1 - 2gd - igH& , V = 1 + ig, a X - 6gH& (58) 

for our model, which can give the sound speed squared c 2 s — C/T>. One can also get the power spectrum and non- 
Gaussianities of curvature perturbations by making use of the explicit form (j54[) in the corresponding formulae that 
has been derived in the above sections. 

For later convenience, let us first introduce some more "slow-variation" parameters, namely 

* s ir-' a "4< p s ^ t s -H- (59) 

Ha Q gti g t<T H g.aatl 

We can always choose the form of g(a) such that the last three parameters \a\, |/3|, I7I 1 all the time. However, 
I77I can only be small when the field remains slow-rolling. This is important because in curvaton scenario, curvature 
perturbations are produced after the end of inflation, and in that case, the slow-rolling of the curvaton field is not 
always satisfied. 

Now let's turn on to the two cases of generating curvature perturbations one by one. In the first case the curvature 
perturbation is generated when curvaton dominates the universe. That requires the decaying of the curvaton is slower 
than that of radiation which is transferred from inflaton. For instance, it is reasonable to assume that the curvaton 
field is still slow-rolling, and exit in a few number of efolds, so it will not lead to another period of rapid acceleration 
26]. Substituting (fM)) into Eqs. ([57]). PT|) and (|52"j) and letting all the slow- variation parameters defined in (f5T)| be 
small, we can get: 



?>P - (7Z?L)(^o) 2 , (60) 



f local 

Jnl 



requil 
JNL 



(I) _ 5 (3agH& -6egH& + e) 
C ~ 6 [ (l-3gH& ) 1 



(61) 



(') ^ WgH&o (3-e)gH& -l 17 13 2 

C ~ ~2^ { l-3gHa a )(1 + 2^ " W<To) ' (62) 



Similarly, in the second case the curvature perturbation is generated when the curvaton decays and becomes 
equilibrium with radiation, and thus the energy density of curvaton has the same scaling as that of radiation with the 
ratio r = p a /pr- Again using (f5"4"]) with Eqs. (fBT)]). (|4"2l and (|53")l , letting a, (3 and 7 be small but retaining r\ for the 
reason given above, we can get: 

^T^(i^p)[( 3 - e + ^)gH& -\{r } + 2,)f{H/a Q )- 2 , (63) 
20, r (-6(3-e + 2^o + (, + 3)) ]W . o)2) (g4) 





flocal 

Jnl 


(II) 


c ~ 


requil 
JNL 


(II) 


c 



(r + 1) L (3(3 - e + 2rj)gH& - (r] + 3)) 2 

±^ [(3 - e + 2 , )9 H, -i„ ^ 



5(r + 1)A [(3 _ t + 2v)gH . o _ 1 ( , + 3))(1 + U _ U ) _ (65) 
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From above we can see that another variable that appears repeatedly is gH&Q. Defining x := gH&o, and in order to 
avoid ghost and gradient instability, x should be smaller than unity. We consider two asymoptotic cases: i) x could 
be positive or negative, with \x\ <C 1 and ii) x is negative, with \x\ ^> 1, which means that the Galileon term plays a 
unimportant /important role at the end of inflation respectively. The spectrum and non-Gaussianities of the curvature 
perturbation can be reduced as: 
i) \x\ « 1: c 2 ~ 1 



V. 



(ii) 



(^l)(H/ao) 2 , AT 



(1 



V A— 2 



(I) 



Hi Ail + if 



144 Hi: 2 ' (H/a ) 2 



(a 



6 

clocal 



I NL 



equil 
NL 



125x 



486 
20 e (g/do) 2 
(r + !)(?? + 3) 



(H/a ) 



requil 
' JNL 



(II) 

c 



125(r + l)x 
5832 



fa + 3) 



(66) 
(67) 



ii) |z| » 1: c 2 



(r? is small for Case I) 



,(/) 



1 



^(-gxff/cvo) 2 , /fcf (J 



(1 + T-) 2 |X| 



. V / |_2 / 



C 6 
e + 2ry) 2 



(2e - a) , 

(-r-T) 



C 



requil 
JNL 



32 y („ + 2)3 V 47T aV (#/<r ) 2 
(") 5(r + 1)(3 - e + 2^(103 + 118*7 + 4?7 2 ) 



rlocal 
) JNL 



11664|a;|(2 + ?7) 2 



C 



35(3 



17496 
40e(H/& ) 2 



(H/& ) 2 



3(r + l)(3 



■2 V )\x\ ' 



(68) 



(69) 



respectively. 

From the above results we can have a couple of comments on the perturbations of our G-Curvaton model. Besides 
the slow variation parameters which are roughly of order 1, now we have three more free parameters, namely the 
value of x, H and ctq &t the end of inflation, to determine and Jnl- Considering the constraints from CMB that 
the amplitude of power spectrum to be nearly 10 -9 , we still have large parameter space to have considerable large 
non-Gaussianities. For instance, for Case I the curvaton field is still slow-rolling at the end of inflation, where \&q\ 
is small compared to H. If it satisfies | <xo j ~ 10~ 3 H, then from the power spectrum we have H ~ 10 -8 , which is 
slightly lower than that of chaotic inflation. For subcase where \x\ <C 1, we can have /J^" ~ 10 2 just by requiring 



1.91 



10 13 . For subcase where \x\ 1, however, we can have / 



equil 
NL 



10 by requiring \g\ 



10 1 



Moreover, in 



both cases //^£ a remains of order unity. For Case II, the energy density of o is comparable to H , which roughly gives 
<r 2 [l + 0(l)x] ~ H 2 . In subcase |x| <C 1, we have |<to| ~ H, which gives H ~ 10 -5 in order to meet the constraint on 
power spectrum. This in turn gives small /Jyx* with /^£ ai roughly of O(10). In subcase |x| 3> 1, however, we obtain 



10 H, it will lead to H ~ 10 from constraint on power spectrum, 



10 7 . Then the nonlinear estimator /^£ ai is roughly of O(10), and AAA 1 



H. If furthcrly it satisfies, i.e., |<to| 
which in turn gives <tq ~ 10 -3 as well as \g\ ■ 

becomes of order unity 4 . Moreover, due to the non-linear effects of the Galileon term G(X, a)Oa in our model, it can 
be expected that large tensor-to-scalar-ratio can be generated as well as large non-Gaussianities [26j, which can be 
distinguishable from the standard curvaton model, of which large non-Gaussianities are also accompanied with lower 
energy scale of inflation, which leads to small tensor-to-scalar ratio (cf. the second paper in [Hj]). 



IV. CONCLUSION AND DISCUSSIONS 



In this paper, we investigated G-Curvaton scenario, where the curvaton field is acted by Galileon action. This 
opens a new access of curvaton scenarios that the fluctuations could be affected by nonlinear terms. After reviewing 
the standard curvaton mechanism, we started from the action of Galileon field and calculated the spectrum of the 
field perturbation. The power spectrum is suppressed by T> given in (|26l) . rather than Kx m the normal single field 
case. We studied the generation of curvature perturbations in both two possible cases, and apart from curvature 
perturbation, we obtained non-Gaussianities of both local and non-local types. 

We have shown, in this work, that there is large possibilities for G-Curvaton to have consistent power spectrum and 
large non-Gaussianities. We presented a concrete model of G-Curvaton as an illustration. With proper choice of the 



4 In our analysis we also made the approximation of H* ~ H , where is the value of Hubble parameter when the fluctuations of a exits 
the horizon during inflation. 
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parameters, the nonlinear estimator could be made of O(10 2 ). However, because of the large parameter space, the 
conclusion is still highly model dependent. We can expect future observational data to have more rigid constraints on 
the G-Curvaton scenario, for example, if future observations can observe large non-local non-Gaussianity compared 
to local one, more or less we can say that it might due to some nonlinear effects such as Galileon. 

Moreover, since the Galileon field can violate the NEC without the ghost and gradient instabilities, G-Curvaton 
can naturally incorporate a model of curvaton with NEC violation, which might be interesting for studying. Here, 
we focus on the inflationary background, however, in principle, G-Curvaton can also be embedded into alternative 
models to inflation, which will make the curvature perturbation induced in corresponding model have more fruitful 
predictions. For phenomenological aspect, there are also a couple of implications that merits thinking of: one example 
is that large non-Gaussianities could be accompanied with large tensor-to-scalar ratio as well, which is different from 
usual curvaton models and thus can be used as a distinguishment observationally; another is that due to the non-linear 
effects from Galileon term, the decaying process of the curvaton might also be modified, which may in turn change 
the speed and amount of the generated products during reheating 29] 5 . We hope that some upcoming works could 
gain more clear insights into G-Curvaton scenario. 
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